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Abstract 

The borocarbides i?Ni2E>2C (R=Gd, Ho, Er) exhibit a large variety of magnetic states and as 
a consequence rich phase diagrams. We have analyzed the nature of these states by specific heat 
investigations. The data were measured down to 0.5 K and up to 80 kOe. The overall evolution 
of each C mag (T, H) curve is observed to reflect faithfully the features of the corresponding H — T 
phase diagram. Within the lower ranges of temperature and fields, the calculations based on 
linearized field-dependent spin-wave theory are found to reproduce satisfactorily the measured 
Cmag(T, H) curves: accordingly, within these ranges, the thermodynamical properties of these 
compounds can be rationalized in terms of only two parameters: the spin-wave energy gap and 
the stiffness coefficient. For the intermediate fields ranges (Hi < H < H sat ) wherein successive 
field- induced metamagnetic modes are stabilized, the evolution of C mag (T, H) is discussed in terms 
of the Maxwell relation (dC mag /dH) T = T (d 2 M/dT 2 ) H . For the particular case of GdNi 2 B 2 C 
wherein the anisotropy is dictated by the classical dipole interaction, C mag (T, H) across the whole 
ordered state is numerically evaluated within the model of Jensen and Rotter [PRB 77 (2008) 
134408]. 

PACS numbers: 
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I. INTRODUCTION 



The H—T magnetic phase diagrams of the heavy members of the intermetallic i?Ni 2 B 2 C 
(R= rare earth) magnets are characterized by a cascade of metamagnetic transformations, 
higher Neel points T/v, and stronger saturation fields, H sat (see Refs.l, 2 and references 
therein). Both T/v and H sat are well scaled by de Gennes factors testifying to the involve- 
ment of the indirect exchange coupling mechanism.- Moreover, except for GdNi 2 B 2 C, there 
are strong anisotropic forces and the competition between these forces and the exchange 
interaction is one of the main driving mechanism behind the characteristic features of their 
magnetic phase diagrams.-^ In general, the role of these competing forces is analyzed in 
terms of a Hamiltonian that consists of bilinear exchange, single-ion crystalline electric field 
(CEF), and dipolar interactions.-^ Calculations based on such a Hamiltonian yield H—T 
phase diagrams that are in reasonable agreement with experiments on, say, HoNi 2 B 2 C,-^ 
GdNi 2 B 2 C, 7,8 , and ErNi 2 B 2 C. & 

The zero-field magnetic structures of the i?Ni 2 B 2 C compounds can be divided into 
two broad classes: 9 (i) the equal-amplitude, collinear, and commensurate AFM structures 
(i?=Pr, Nd, Dy, Ho ): here the moments are coupled ferromagnetically within the layers 
which are, in turn, antiferromagnetically stacked on each other along the c axis and (ii) the 
amplitude- modulated, incommensurate structures (R= Gd , Tb, Er, Tm); on lowering the 
temperature, these modulated structures are gradually transformed into an equal-amplitude, 
squared-up spin-density states. 

On applying an external magnetic field along the easy-axis, both the commensurate 
and incommensurate magnetic modes are observed to undergo a cascade of magnetic phase 
transformations leading to a set of different magnetic modes; each having a distinct prop- 
agation vector and a distinct total magnetic moment. For the case of, say, HoNi 2 B 2 C and 
ErNi 2 B 2 C, there emerges a very rich and interesting H-T magnetic phase diagram^^^ 
with characteristic features that are reflected in all measured thermodynamical proper- 
ties. As an example, the low-temperature magnetization isotherms show a cascade of step- 
like metamagnetic phase transformations with constant plateau between the consecutive 
transformations.— iMyULiLIZ As the low-temperature magnetic properties are governed by the 
excitation spectra of the magnon quasiparticles then, according to the spin-wave theory, 
the surge of such a stair-like features is a manifestation of related field-dependent distinct 
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features in the dispersion relations or density of states of the magnon quasiparticles. Then a 
detailed spin-wave analysis of the H- and T-dependence of the thermodynamical properties 
would contribute positively to the understanding of the magnetic phase diagrams and the 
related thermodynamical properties of these borocarbide magnets. Using field-dependent 
specific heat measurements on single-crystals of i?Ni 2 B 2 C (R = Gd, Ho, Er), this work in- 
vestigated the magnetic contribution to the specific heat when T and H are varied across 
these H-T phase diagrams. 

The magnetic structures of the compounds under study are well established for the high- 
and low-field part of the phase diagram: within the intermediate-field range, the structures 
are either unknown or too complex. Furthermore, it happened that the highest magnetic 
field available for this study is lower than the one needed to reach the high-field saturated 
range. Accordingly, the comparison of the measured and calculated magnetic specific heats 
(the latter is based on spin-wave theory) was carried out only for the lower field limit (§ 
II. A and IV. C). For the intermediate-field range (as well as across the available H — T plane 
of the phase diagram), basic thermodynamical relations were employed to relate the field 
evolution of the measured magnetic specific heat to the thermal evolution of the isofield 
magnetization curves (§ II. C and § IV. C). Finally, since the linearized magnon theory is 
applicable only for the lower temperature range, mean-field model calculations were carried 
out to evaluate the magnetic specific heat within the higher range of the ordered state; this 
part of the study was restricted to only GdNi 2 B 2 C (§ II. B and § IV.C.l.b): in contrast to R 
= Ho and Er, GdNi 2 B 2 C has negligible crystalline field anisotropy and as such the classical 
dipole interactions play a crucial role in shaping its magnetic structure. 

II. THEORETICAL BACKGROUND 

A. Field-dependent magnon specific heat of iiN^E^C. 

The magnetic properties of i?Ni 2 B 2 C are usually described in terms of the Hamiltonian: 



where Tiz is the Zeeman term with the field H applied at angle 6 away from the easy axis: 



T~L = T~tz + 'Hex + 7~Lcef-> 



(1) 




(2) 
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where Sj is the total angular momentum operator of the j th R 3+ -ion, gj is the Lande factor, 
and He is the Bohr magneton. The bilinear isotropic exchange interaction is: 



Hex = -\Y,Jft)%£j> ( 3 ) 

where is the isotropic coupling between Si and Sj. Considering the low-temperature 

commensurate AFM structures of i?Ni 2 B 2 C then 7i ex (Eq. [3]) can be split into two terms: 

Hex = + 22 ^ij^Sf + ^2 ^ij B St-Sf ( 4 ) 

<ij>,LeA,B <ij>A,B 

j£ couples moments S[ and Sj' within the same layer (L = A or B) and J^ B couples 
moments from different layers. At lower temperatures, we assume that the net effect of H.cef 
on the magnetic moment can be approximated by an easy-axis anisotropic field (represented 
by H a ),—^2- consequently: 

Hcef + n z = -gnB(H a +H).J2sf~ 9VB(H a - H). ^ Sf . (5) 

igA jeB 

Let us consider the following useful simplification: as far as the magnon propagation is 
concerned, the low-temperature squared-up incommensurate magnetic structures is assumed 
to behave as if it is a collinear AFM structure (for, say, ErNi 2 B 2 C, this ignores the presence of 
the kinks along the a axis): such a simplification is motivated by the remarkable feature that 
though i?Ni 2 B 2 C ( i?=Ho, Er) have different magnetic structures, - however they manifest 
similar stair-like features in their M(T, H) isotherms^ 1 ^^^ as well as similar expression 
for their zero- field C mag (T, H) curves.— 

Using standard linearized spin-wave procedures,— the total Hamiltonian (Eq. [1]) is diag- 
onalized and consequently the magnetic specific heat is calculated. Based on the strength 
of the applied field, three different situations can be distinguished: (i) the low-field limit 
(H < H s f. the spin-flop field) wherein the spins are basically aligned along the preferred 
orientation, (ii) the intermediate field regime wherein a series of field-induced metamagnetic 
states are stabilized, and (iii)the high field limit (H ^ H sat : saturated field) wherein spins 
are driven towards paramagnetic saturation. 
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1. Low-field approximation (H < H s f) 



The dispersion relation of the two modes (+ and -) are: 



tvjj± = +VA 2 - B 2 ± g/2 B Hcos(e), 

A = AS Jo {1-0.5 [cos{ak x ) + cos(ak y ))} + 85 Ji + gn B H a , 

B = 2SJi |cos(-fc x + -k y + -k g ) + cos(-fc x + -k y - -k z ) 

,a a c a a c . i 

+ cos(-A; x - -k y + ~k z ) + cos{--k x + ~k y + -k z )j (6) 

where j7o and J7i are effective exchange couplings among the nearest neighboring within, 
respectively, the same plane and within neighboring planes. Accordingly, there are two 
distinct field-dependent energy gaps: 



A±(#) = +^( 9 fx B H a ) 2 + 16 J l9 fx B SH a ± 9 fi B Hcos(9) = A a ± gfx B Hcos(9); (7) 

evidently A ± (H) do not depend on J and their separation is linearly related to Hcos{9). 

At lower temperature where long-wave limit is valid, the integration over the allowed 
fc-space yields the magnon specific heat: 

oo 

C mag (T, H) = RA 4 J(47r 2 D 3 a T 2 ) ^ {T cosh(m£/T) [K 4 (mA a /T) + <1 + 2? / A 2 a ) K 2 (mA a /T 

m=l 

-0.5m£ sinh(m£/T) [K 4 (mA a /T) - K (mA a /T)}} (8) 

where K n (mA a /T) represents the modified Bessel function of the second kind, £ = 
gix B Hcos{6) /k B and D a is a measure of the stiffness and is function of the exchange couplings 
and magnetic anisotropy: 

D a = (16(Jq + Ji)J S 2 + 2J S( 9f i B H a ))K (4 J X S)* (9) 

Due to the type of the undertaken approximations, Eq. [8] does not hold for H cos(6) 
equal or higher than the spin-flop field H s f which is the value at which the lower branch 
goes to zero. The field influence on Eq. [8] enters only through Hcos{9). Furthermore, this 
equation reduces to Eq. 25 of Joenk 19 when Jq = J\ and it gives the well-known T 3 relation 
when k B T > A a .i£ 

For T — > 0, A a /T — > oo, and £ < A , Eq. M reduces to: 



RAl /2 X-A a , 



CmagiT, H) ~ 21/2^3/2 eXp ( J> 

indicating a dominant exponential character within this H — T region. 



A n \A n 



(10) 



2. Intermediate fields range (H s f < H < H sat ) 

As mentioned in § I, the absence of a detailed description of the involved magnetic 
structures together with the absence of analytical expression for the dispersion relations 
hinder any direct evaluation of the magnon specific heat contribution within this range. 

3. The high- field limit (H > H sa t) 

Due to the limitation of our experimental conditions, we were not able to probe the 
magnon contribution within the region of the H — T phase diagram wherein the induced FM 
state is established. Nevertheless for completeness sake we derive the magnetic contribution 
for this state in appendix A below. 

B. Mean-field Model calculation of C mag (T, H) of GdNi 2 B 2 C 

The above calculations of the magnon C mag (T, H) of i?Ni 2 B 2 C compounds are valid only 
within the lower temperature range, mostly below liquid helium temperatures. To inves- 
tigate the magnetic contribution within the whole magnetically ordered range, we resort 
to mean-field model calculations applied on the simpler case of the ordered 7 ' 2 S -moments 
of GdNi 2 B 2 C. Jensen and Rotter- showed that a model consisting of a sum of the bilinear 
Heisenberg exchange term and the classical dipole interaction is able to explain the features of 
the whole magnetic phase diagram as well as the so-called magnetoelastic paradox.— and the 
zero-field C mag (T) of GdNi 2 B 2 C. here, in this work, their calculations of C mag (T, H) are ex- 
tended to fields up to H < 80 kOe using the McPhase program package (www.mcphase.de).— 

C. Application of the Maxwell relation to correlate M(T,H) and C mag (T, H) 

An applied field on an AFM mode tends to remove the degeneracy appearing in Eqs. 
[6] and UJ As a consequence, the fractional contribution of the lower mode to C mag (T, H) 
increases leading to (dC mag /dH) T > 0. On the other hand, for H — > H sat , an increase in 
H would induce a gradual decrease in C mag (T, H) resulting in (dC mag /dH) T < 0. Similar 
arguments hold for M(T,H). In fact the following Maxwell equations relate the evolution 
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of M{T, H) to that of C mag (T, H): 



(dS mag /dH) T = (dM/dT) H 
(dC mag /dH) T = T (d 2 M/dT 2 ) 



(lib) 



(11a) 



These general relations are very helpful in the analysis of S mag (T, H) and C mag (T, H) since 
it is much easier to measure the thermal evolution of M(T, H) than the field dependence 
of C mag (T, H). Furthermore, they hold across the whole H — T phase diagram and are 
independent of field orientation (whether H\\a or H\\c) or the type of the involved magnetic 
structure (whether commensurate or incommensurate). This utility is most welcomed when 
investigating those regions of the H — T phase diagrams wherein the spin- wave analysis of 
§ II. A is not applicable. 

Evidently each of C mag (T < T N , H) and M(T < T N , H) would mirror the exotic features 
of the H — T phase diagrams of i?Ni 2 B 2 C;— iMii£ii£ii& in fact, these diagrams have been 
constructed from the magnetic anomalies occurring in C mag (T < T N , H) and M(T < T N , H) 
as well as many other magnetic properties such as the integrated intensities of the magnetic 
neutron diffractograms;— ^ the magnetostriction,— 1 ^ and the magnetoresistivity.— ^ Here 
in this work, we use Eqs. Illal and lllbl to discuss the general trend of C mag (T, H) in terms 
of the trend of the extensively reported M(T, H) curves.- 

III. EXPERIMENTAL 

The i?Ni 2 B 2 C (R— Gd, Ho, Er) single-crystals were selected for this study because 
they offer a good representation of the magnetic properties of the whole i?Ni 2 B 2 C series: 
HoNi 2 B 2 C is a typical representative of the collinear, commensurate, AFM structures while 
ErNi 2 B 2 C and GdNi 2 B 2 C are good representatives of the modulated, incommensurate mag- 
netic structures; the anisotropic forces are stronger in the former— while extremely weaker 
in the latter.— GdNi 2 B 2 C has negligible CEF forces; nonetheless, it is shown that the 
dipolar forces are essential for directing the moments transversal to the propagation vector 



The single crystals of these three representatives, together with that of the reference 
YNi 2 B 2 C, were grown by floating zone method.— Results from extensive structural and 
physical characterizations are in good agreement with the reported data, confirming the 
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good quality of our crystals. The temperature-dependent specific heat at fixed fields was 
measured on two different setups. One is pulse-type adiabatic calorimeter [500 mK <T<25 
K, f 20 kOe] and the second is a quasi-adiabatic setup with a temperature range covering 
f .5-1 00K and a field up to 80 kOe. For the sake of completeness, we quote some of our results 
(in particular the zero-field magnetic specific heat) that had been reported in Refs.— i22i2Q 

IV. RESULTS AND DISCUSSION 

A. Electronic and lattice contribution 

For all compounds, the total specific heat C to t was analyzed as a sum of an electronic C e , 
a phonon C p h, a nuclear Cjy, and a magnetic contribution C mag : 

Ctot = Ce + Cph + C n + C mag (12) 

For consistency reasons, the lattice contribution was estimated from the specific heat of the 
isomorphous single-crystal of YNI2B2C.— We tried other means of estimating the phonon 
contribution (such as mass normalization 3 -^ or using other nonmagnetic isomorphs 29,32 ). It 
is found out that the spin-wave fit-parameters based on different methods are differing by 
only a few percent; other than this variation, the use of different estimation process does 
not influence the conclusions drawn from this work. 

The low-temperature, normal-state electronic contribution of these compounds was eval- 
uated from C to t(T, H) of YNi 2 B 2 C single-crystal (the same as that used in Ref .— ). Within 
the normal state (or H > if applicable), the electronic, phonic, or nuclear contributions 
are taken to be field independent (see below). An illustration of the various contributions 
of i?Ni 2 B 2 C (R =Ho, Er) are shown in Fig. [TJ evidently for all R— Gd, Ho, Er compounds, 
the diamagnetic contribution is small in comparison with the magnetic or nuclear term. 

B. Nuclear Contribution 

The magnetic contribution was obtained as follows: after subtracting the electronic and 
phonon by the process explained in §IV.A, the resultant {C to t — C e — C p h) is confronted 
with the sum of the nuclear (Eq. [T3l) and magnetic (Eq. [8]) terms. Finally, the magnetic 
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FIG. 1: (Color online) Log- log plot of the zero- field Ctot{T) versus T showing the various individual 
contributions that are contained in Eq. [T2l (a) H0M2B2C and (b) ErNi2E>2C. The triangles (cir- 
cles) denote the experimental total (magnetic) specific heats, while the lines represent the various 
calculated contributions: dotted, dashed, solid-thin (red) and solid-thick (black) lines represent, 
respectively, the electron+phonon, nuclear, magnetic, and total contribution. Both nuclear and 
magnetic contributions are obtained from the least-square fits (see text). 

contribution is obtained after subtracting out the nuclear specific heat term which is usually 
given as:— ^ 



/ +/ +1 



— L)j — L)i 



) / EE ex p( 



\ r +/ +/ 



— UOj — UJi 



C n {T) = A tso J2 E K 



) exp( 



k B T 



) 



\i=-ij=-i 



J li=-ij=-i 




(13) 
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where A iso is the isotope abundance, a int and P are the magnetic dipole and electric 
quadrupole interaction parameters of the nuclear spins, respectively; I is the total nuclear 
spin while i is its component along the quantization axis. a ex t is associated with the exter- 
nally applied magnetic field which is extremely small if compared to the internal field: for, 
say, ErNi 2 B 2 C,— H int ~7xl0 6 Oe and thus the highest applied field is only ~1% of H int . 

TABLE I: The nuclear hyperfine parameters (Aj so , cti n t, and P) and the spin- wave parameters 
(A a and D a ) of H0M2B2C and ErNi2B2C which are obtained after fitting the experimental 
(Ctot — C e — C p h) curves to the sum of the nuclear (Eq. [T3j) and magnetic (Eq. [8]) terms. The 
specified temperature range indicates the region wherein the nuclear contribution is dominant (for 
the temperature range of the fit see § IV. C. 2 and § IVC.3). The obtained nuclear parameters are 
compared with the corresponding parameters of the R- metal (see text). 



i?Ni 2 B 2 C 


Temp. ran£ 


;e A iso I aint PA D 




(K) 


% K mK K K 


HoNi 2 B 2 C 


[0.5,1.5] 


100 7/2 0.362 -9.6 7.7 ±0.3 4.6±0. 2 


Ho metal 33 


[0.03,0.5] 


100 7/2 0.32 7.0 


ErNi 2 B 2 C 


[0.1,0.5] 


90(5) 7/2 0.054 2.9 7.0 ±0.1 3.0±0.1 


Er metat" 


[0.03,0.8] 


100 7/2 0.042 -2.7 



The least square fit involves the simultaneous search for the best values of the five param- 
eters, namely the nuclear parameters A iso , ai nt , and P as well as the spin- wave parameters 
A a and D a (see below). After the substitution of the obtained nuclear parameters (Table 
HJ, the expression of C n (T) (see Eq. [13]) does reproduce satisfactorily the measured nuclear 
specific heat of ErNi 2 B 2 C and HoN^E^C (GdN^E^C has no nuclear contribution): the 
overall fits are shown in Fig. [IJ It is evident that a int of the studied R 3+ compounds are 
close to the values reported for the corresponding rare-earth metala 3 ^^ indicating that the 
hyperfine field is determined mainly by the internal electronic configuration of the R 3+ ion. 
On the other hand, the P parameters are extremely small; this is not surprising since the 
point group of the sites at which the R 3+ nucleus resides is in borocarbides and in 
the elemental rare earth. 
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C. Magnetic contribution 

Based on the general features of the C mag (T, H) curves (shown in Figs. [21 El and ED, 
one distinguishes four temperature regions: (i) a paramagnetic region, T > Tn(H), wherein 
C m ag{T, H) is due to change in the population of the crystal field levels, (ii) a critical region, 
T m Tjf(H), wherein C mag (T, H) is related to critical phenomena, (iii) an intermediate 
region, T x < T < T N (T x = T R for GdNi 2 B 2 C or T WFM for ErNi 2 B 2 C) which encompasses 
the sine modulated states. Within this region, the spin- wave analysis of § II. A. 1-2 is not 
applicable. Finally (iv) the low-temperature region which should be restricted if one intends 
to analyze C mag (T, Hp) in terms of the linearized magnon theory: below we restrict this 
range to T < A and H < H s f (none of the compounds under study presents the H > H sat 
case). Within this region, the collinear AFM/squared-up states are established and the 
measured C mag (T, H) is to be confronted with Eq. [HI It is worth mentioning that for 
extremely low-temperatures (T << A) , the exponential decaying character of C mag (T, H) 
(see Eq. [TOl) and the relatively large contribution of the nuclear Schottky contribution limit 
the usefulness of the lower temperature range for the magnon analysis. 

1. GdNi 2 B 2 C 

The zero-field magnetic structure of GdNi 2 B 2 C (T x = 19.5 K) is an incommensurate 
sine-modulated structure (moments along b axis and ~q =0.551a*).— At Tr ^13.5 K, a 
moment reorientation sets-in leading to an additional modulated mode transversely polar- 
ized and having a small amplitude along the c-axis.— ^ Two H — T phase diagrams were 
reported:— one for H || a and another for H || c. The phase diagram for H\\a shows three 
field-induced magnetic phase transitions [see inset of Figj2] (a)]: (i) the saturation boundary 
with H sat ; (ii) the reorientation boundary H R a ; and finally (iii) the domain- wall boundary 
H^. In contrast, the phase diagram for H\\c shows only two transitions [see inset of FigfS] 
(b)]: HgQ t and H R C ; in comparison, the thermal evolution of the the former (latter) is similar 
to (different from) that of i?J" t (H R a ). None of the two phase diagrams shows those char- 
acteristic i7-induced cascade of metamagnetic phase transitions which are common in, say, 
the case of R = Er, Ho; this is attributed to the absence of strong anisotropic features. 

As mentioned above, the experimental magnetic specific heat of GdNi 2 B 2 C would be 
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confronted with two model calculations: (i) the magnon calculation which is valid for sub- 
helium temperature (§ II. A) and the model calculations of Jensen and Rotter- which are 
extended from the helium-temperature range up to Tjy (§ II. B). 



FIG. 2: (Color online) The isofield C mag (T, H) of GdN^E^C at various magnetic fields applied 
along (a) H\\a and (b) H\\c. The inset at the bottom-left (top- left) illustrates the H — T phase 
diagram for H\\a (i/||c); 23 there, the horizontal dotted lines represent the fields that were applied 
during these measurements. Evidently, both and anomalies in the C mag (T, H) curve do 
reproduce the boundaries of the phase diagrams, in particular the reentrant feature appearing for 
the H\\a case (see text). 

a. Magnon contribution to C mag (T, H) Figure [21 shows that the most prominent fea- 
tures of C mag (T, H) of GdNi 2 B 2 C are (i) the characteristic and distinct evolution of the 
HgatiT) and H R (T) curves and that (ii) within a certain region of T and H, C mag (T, H) 
appears to be hardly influence by H. The cause of this apparent collapse of the C mag (T, H) 
curves becomes clear if we compare these curves with the predictions of Eq. [HJ Fig. [3] 
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shows a fit of the measured C mag (T, H = 0) curve to Eq. [H]and the obtained fit parameters 
are A = 2.9 ±0.1 K and D =5.6±0.1 K which are close but, due to the difference in the 
temperature range, are better than those reported in Ref.— . As both A a and D a are field 
independent, then the insertion of these values into Eq. E] leads to the calculated C mag (T, H) 
curves (with no adjustable parameters) for all fields up to H < ksA a / (g^is) ~ 20 kOe which 
is the field above or equal to which the magnon calculations based on Eq. [S] are not valid. 
Evidently the calculated and measured curves collapse on each other when T/A a is in the 
immediate neighborhood of 1: thus the apparent collapse is a reminder that our experimen- 
tal conditions are good only for probing that part of the phase diagram wherein the field has 
a weak influence on the strongly exchanged-coupled AFM-like state. It is worth mentioning 
that such a collapsing feature is reflected also, by virtue of Eq. 111b} in the low-temperature 
M(T, H) curves^ 

Figure [5] shows also that the if— evolution of the low-temperature C mag (T, H\\ c ) curves is 
similar to that of C ma g{T, Hu a ) ones; the only difference is that all values of H\\ a are lower 
than H R a {T < Tr) boundary while, in contrast, some of the applied H\\ c values are higher 
than H R C (T < Tr) (see above); then it is no surprise that the presence of the Tft(f7|i c )-event 
in C mag (T, H\\ c ) is more pronounced than in C mag (T, Hn a )\ in fact there is no manifestation 
of the reorientation event in the isofield C mag (T, H\\ c >32kOe) curves while, in contrast, 
for H\\a, the presence of the reorientation event is evident in all applied field up to the 
maximum 80 kOe (see insets of Fig. [2]). It is noted that the reentrant feature of the H R a (T) 
curve within the neighborhood of 11 K is well evident in the C mag (T, H\\ a ) curves: as H\\a is 
increased, the peak associated with this event moves first to higher temperature but reverts 
to a decreasing tendency when Hu a reaches values higher than 30 kOe. 

Within T < 5 K and ff„ a < H R (T) range, (dS mag /dH) T , (dC mag /dH) T , and (dM/dT) H 
[Ref.-22a2Z] are weak but positive. Similar features are evident for the H\\c case. When 
H\\ c — > Hr(T), C mag (T, H\\ c ) is observed to increase, reaching a maximum at the phase 
boundary. A further increase in H\\ c > Hr(T < Tr) leads to {dC mag / dH)r < indicating 
a decrease in the entropy and as such an increase in the ordered component along the c 
axis. For T > T R , both (dM/dT) H [Refe/M] and {dC mag /dH) T are weak and positive for 
H <30 kOe but negative for 30 kOe< H < H sat . 

b. Model calculation of C mag (T, H) Figure H] compares the measured magnetic specific 
heat of GdNi 2 B 2 C with the model calculation of Jensen and Rotter (see § II. B).- It is assuring 
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FIG. 3: (Color online) Log-Log plot ofC mag versus the normalized temperature (T/A) under differ- 
ent applied fields. A a represent the zero-field energy gap of Eq. [7J Symbols denote measurements 
while solid lines represent the calculated C mag (T, H) based on Eq. [8j The zero-field C mag (T, H) 
curve (measured down to 0.5 K) was fitted with Eq. [HI The obtained parameters were fed into Eq. 
[8j and there from the theoretical C mag (T, H) curves for field up to 20 kOe were calculated with no 
adjustable parameters (see text). 

to notice that, even though there are no adjustable parameters in these calculations, the 
model is able to reproduce the main features of the measured C mag {2 K< T < T N ) for both 
H\\ a < 80 kOe and H\\ c < 80 kOe [see, respectively, Figs. H^a) and (b)]. The following three 
achievements of the model calculations should be highlighted: (i) the calculated magnitude 
of the steps at both Tr and T/v compare favorably with the measured values; (ii) the surge 
of an anisotropy for the spherical 7 I 2 S Gd- moments even at temperatures as high as 20 K is 
well accounted for; and (iii) the reorientation process at Tr (along both field orientations) is 
shown to be a consequence of the joint action of exchange interactions and dipolar forces even 
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FIG. 4: Comparison of the calculated (solid lines) and measured (symbol) field-dependent magnetic 
specific heat of GdN^E^C for (a) H\\a axis and (b) H\\c axis. For ease of visualization, the 
successive curves are displaced upwards with the same amount of vertical shift (10 J/moleK). 
The calculations (with no adjustable parameters) are based on the model proposed by Jensen and 
Rotter— (see text). 

though the energy of the former is at least five times larger than that of the latter: isotropic 
bilinear interaction, alone by themselves, do not lead to any reorientation processes.— i2L2£ 
Evidently in spite of the above-mentioned successes, this mean-field model calculation are 
not expected to account for the magnetic features of the specific heat within the very low- 
temperatures or in the neighborhood of T N (the critical region). 
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2. HoNi 2 B 2 C 



The zero-field magnetic ground structure of HoNi 2 B 2 C is an AFM state wherein the 
ferromagnetic basal planes are piled along the c axis forming an AFM configuration /* / 
with q = (0,0,1).— On applying a field along the easy (1,1,0) axis, the phase diagram 
of HoNi 2 B 2 C shows a succession of field-induced metamagnetic phases,— transforming 
from '/* into: at Hi, /* y\ at H 2 , and a saturated /'/'/' state at H sat . On 

the other hand, if the field is applied within the ab plane but at an angle 9 away from the 
easy axis, then the strong CEF-induced anisotropic character of HoNi 2 B 2 C becomes evident 
and as a result restricts the moments to be along only the <110> directions: as shown 
in the lower inset of Fig. [5j an increase in 9 induces an increase in both Hi [expressed as 
4.1/cos(0) kOe] and H sat [related by 6.6/sin(45-6 l ) kOe] but a decrease in H 2 [described by 
8.4/cos(45-fl) kOe]M 

Figure shows the thermal evolution of the isofield C mag (T, H) of HoNi 2 B 2 C for various 
fields, all applied along the (100) axis (9=45°); its overall field-dependence can be compre- 
hended in the light of Eq. EB the M(T, H) curves, and the features of the magnetic phase 
diagrams. One of the most remarkable features of these C mag (T, H) curves is the observation 
that their field-dependence (the upper inset of Fig. Ej) is reminiscent of the field-dependence 
of the isothermal M(T,H): as an example, for H < Hi,— (d 2 M/d 2 T) H > and concomi- 
tantly C mag (T, H) (see Fig. ED increases with H. Furthermore, C mag (T, Hi < H < H 2 ) is 
relatively much higher but C mag (T,H > H 2 ) decreases monotonically with the field as it it 
approaches the forced saturated /'/'/' state. 

Based on the H — 9 phase diagram (lower inset of Fig. [5]), it becomes clear that among the 
various C mag {T, H) curves, the only candidates to be contrasted with the spin wave analysis 
of § II. A. 1 are those measured within the H <Hi range, namely Hn a = and 2.5 kOe wherein 
the low-T collinear AFM /* y state is established.— Fig. [U] shows the satisfactorily fit of 
C mag (T, H = 0) and the best fit parameters are: A = 7.7 ±0.3 K and D =4.6±0.2 K. These 
values are in agreement with our earlier results.— Based on Eq. [HJ the thermal evolution 
of the C mag (T, if|| a = 2.5 kOe) should be obtainable form these A and D and the value of 
the field component along the (110) direction; this is indeed the case as can be convincingly 
observed in Fig. [HI 



17 



FIG. 5: (Color online) Thermal evolution of the magnetic specific heats of H0M2B2C under various 
applied fields. 30 For clarity, the C mag (T,H < 10 kOe) curves are plotted in the lower panel while 
the C mag (T, H > 10 kOe) curves are plotted in the upper panel. The lower inset shows the angular 
dependence of the critical fields (adapted from Ref.—); the horizontal dashed lines represent the 
applied fields. The upper inset shows the field-dependence of the magnetic specific heat at 2 K: 
C mag (T, H) within the same metamagnetic state is almost the same; this feature is reflected also 
in the M(2K, H) isotherms reported in Ref.— (see text). 

3. ErNi 2 B 2 C 

The magnetic phase diagram of ErNi 2 B 2 C is shown in the inset of Fig. [7J— in zero field, 
this compound superconducts at 10.5 K and orders magnetically at TV =6.4 K into an 
incommensurate modulated AFM structure with g=(0.55,0,0) and moments pointing along 
the b axis.— Below Twfm =2.2 K, a weak ferromagnetic (WFM) stated emerges together 
with an equal-amplitude, squared-up state.— A series of field-induced metamagnetic trans- 
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FIG. 6: (Color online) A log-log plot of the isothermal C mag (T/A a ,H ) curves of H0N12B2C. A a 
represent the zero-field energy gap of Eq. [7J For clarity, the C mag (T/ A a , H < 10 kOe) curves are 
plotted in the lower panel while those of C mag (T/A a , H > 10 kOe) curves are plotted in the upper 
panel. The solid line represents the non-linear least square fit of the zero-field measured curve to 
Eq. E] (see text). 

formations appears when a field (H<H sat ) is applied along the easy axis (0,1,0):— J&ALiSi a t 
2 K, three metamagnetic transformations occur at 4, 11, and 20 kOe. For H > H sat , the 
paramagnetic saturated state is stabilized. 

To investigate the H- and T-dependence of C mag {T, H) across such a phase diagram, 
we carried out a series of isofield measurements with if||(010). The resulting C mag (T, H\\ a ) 
curves are shown in Figs. [71 The thermal evolution of the zero-field entropy (not shown) 
suggests that the lowest four level are fully populated above 10 K; this result is in agreement 
with the findings of Gasser et alM that the electronic ground state is a doublet which is 
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FIG. 7: (Color online) Thermal evolution of C mag (T, H) of ErN^E^C at various magnetic fields. 
The inset shows the H — T phase diagram (adapted from Ref.™): the thick dashed line represent 
the H C 2(T) —superconductivity— curve while the solid lines represent the various magnetic trans- 
formations. The horizontal dotted lines represent the magnetic fields that were used during this 
study. 

separated from the immediate excited state (also a doublet) by 0.6 ~ 0.7 meV (7~8 K). 
For H < Hi ( T < Twfm), an increase in H induces an increase in S mag (H) and, based on 
Eq. Illal and the experimental results,— {DM / 8T)h > (considerations should be made for 
the presence of the superconductivity below H c2 and the surge of weak ferromagnetic state 
below T WFM ). 

For T < T WFM (H < Hx), the zero-field magnetic sate is approximated as an AFM 
structure (see § II. A ); accordingly, the measured C mag (T < Twfm,H < Hi) curves are 
confronted with Eq. [HJ Fig. [8] shows the excellent fit of C mag (T, H = kOe) which gives 
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FIG. 8: (Color online) A log-log plot of the Cm(T,H) curves of ErN^E^C at various applied 
magnetic fields. The magnetic specific heat (symbol) are compared with the theoretical calculation 
(lines) based on Eq. [8] (see text). The calculated curves of H=0, 5, 10 kOe are different from each 
other only at T < A. 

A = 7.0 ±0.1 K and D =3.0±0.1 K. Using these parameters, we calculated the field- 
dependent C mag (T, H) curves (see Fig. [SJ. Once more, one observes the collapse of the 
Cmag(T < T WFM ,H < Hi) curves within the immediate neighborhood of T^ FM . Due to 
experimental limitations, we were not able to probe the field evolution for T < Twfm nor 
the thermal evolution for H > H sat of C mag (T, H) of ErNi 2 B 2 C. 
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V. DISCUSSION AND SUMMARY 



The general evolution of C mag (T, H) within the magnetically ordered state of represen- 
tative .RM2B2C (R=Gd, Ho, Er) is found to reflect faithfully the characteristic features of 
their H — T phase diagrams. Three approaches are employed for the analysis of the evolution 
of these C mag (T, H) curves: (i) basic thermodynamical analysis which allows us to relate 
the evolution of C mag (T, H) to that of the magnetization measurements: thus permitting a 
generalization to field and temperature ranges beyond the limitation of our experimental set 
ups. (ii) the linearized spin-wave analysis which allows us to investigate the low-temperature, 
low-field range; and (ii) the model calculation based on which we are able to probe the higher 
temperature region, a region which is not accessible for spin-wave analysis. 

One of the characteristic features of the magnetic phase diagrams is the stair-like behavior 
observed in the magnetization isotherms: it is assuring that this feature is manifested also 
in the C mag (T, H) case: all the C mag (T, H) curves within the same metamagnetic mode do 
collapse on each other (see e.g. Fig. [5]). This feature is related to the influence of the 
elementary magnetic excitations on the thermodynamical properties and here, in this work, 
we discuss this influence in terms of the linearized spin-wave model (see the theory in § II. A 
and its confrontation with the measured C mag (T, H) curves in § IV.C) 

Based on the analysis of § IV one is able to delineate the low-field, low-temperature range 
of the H — T phase diagrams wherein the linearized magnon contribution is found to describe 
satisfactory the experimental results: for GdNi2B2C, it is the H < Hr(T) range wherein the 
sine- modulated state is squared- up; for ErNi 2 B 2 C it encompasses the H < Hi(T < Twfm) 
range wherein the squared-up AFM-like state is established; and finally for HoNi2B2C, it 
is the H < Hi (T < Tjy) range wherein the collinear AFM structure is established. In all 
these magnetic states, the dispersion relation is taken to be given by Eq. [U]and the magnetic 
specific heat is expressed by Eq. [HJ the satisfactorily agreements between the calculated and 
measured C mag (T, H) justifies the assumptions considered in this model. 

The values of the fit parameters evolves reasonably well across the studied compounds: 
as an example, the gap parameter of GdNi2B2C (A = 2.9 K) is much smaller than the 
corresponding values of HoNi 2 B 2 C (A = 7.7 K) and ErNi 2 B 2 C (A = 7.0 K); based on 
Eq. [71 such a result does agree with the well-established fact that the anisotropic field 
for GdNi2B2C is extremely small. On the other hand, the stiffness constant of GdNi2B2C 
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(D = 5.5 K) is greater than that of HoNi 2 B 2 C (D = 4.6 K) and ErNi 2 B 2 C (D = 3.0 K); 
based on Eq. EH this is related to the fact that the effective exchange couplings of GdNi 2 B 2 C 
(proportional to the de Gennes factor) are the strongest. It is recalled that no direct scaling 
with the de Gennes factor should be expected since the D term of Eq. [9] contains also H a . 

As the investigated compounds are good representatives of the other magnetic borocar- 
bides, then it is expected that the above mentioned characteristic step-like behavior should 
be manifested also in other magnetic i?Ni 2 B 2 C compounds: indeed, the M(T, H) isotherms 
of TbNi 2 B 2 C (having similar anisotropic features as those of ErNi 2 B 2 C) manifest such a 
stair-like feature.— Mdk Then, based on Eq. Illbl its C mag (T, H) features should be similar 
to the those of the studied compounds. It is emphasized that the manifestation of such 
stair-like features is a more general property since it is manifested, not only in the M(T, H) 
isotherms, but also in several other thermodynamical quantities (see § II. B). This generality 
suggests that, although other models have been applied to analyze these stair-like features 
in, e.g., HoNi 2 B 2 C [Ref.-^] and ErNi 2 B 2 C,- an analysis in terms of the spin- wave model is 
shown to be extremely useful for the description of the low-temperature thermodynamical 
properties of these magnetic compounds. 
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APPENDIX A: THE HIGH-FIELD LIMIT (H > H sat ) 

For the induced FM state, we assume that the action of both the anisotropic field and 
exchange couplings is the similar to the case of the low-field limit (§ II. A) except that here 
the field is strong enough to overturn the antiparallel spins: all the spins are oriented along 
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the easy-axis directions; e.g. for the case of HoNi 2 B 2 C, it is one of the four <110> directions 
that is nearest to the field. In the long wave limit, the dispersion relation simplifies to: 



ftw k = A f + Cl (k 2 x + k 2 y )-c 2 k 2 z , 



(Al) 



where, assuming a weaker orthorhombic distortion (a^b): 



ci = 2Sa 2 (J - J x ) 
c 2 = 2Sc 2 ^Ti, 



and the energy gap {k = 0): 



A f = 2gfi B Hcos(9)/k B . 



(A2) 



Thus, according to the above assumptions — involving the anisotropic field and exchange 
couplings — the spin-wave parameters of the saturated FM state are different from the 
standard (spontaneous) FM state: in the latter case the gap includes a 2H a term and 
the sign of the c 2 term in Eq. IA1I is positive. Applying standard procedure, the magnon 
contribution is given as: 



For T < A f , Eq. \M\ reduces to RAexp(-A/T)(ii*l 2 D 3 f /2 T 1 / 2 ) while for T » A f it 
gives the well-known T 3 / 2 relation.— 




(A3) 



where 



D f = 2S{J -J 1 )sj°. 



(A4) 
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